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We study AKLT models on locally tree-like lattices of fixed connectivity and find that they exhibit 
a variety of ground states depending upon the spin, coordination and global (graph) topology. We 
find a) quantum paramagnetic or valence bond solid ground states, b) critical and ordered Neel 
states on bipartite infinite Cayley trees and c) critical and ordered quantum vector spin glass states 
on random graphs of fixed connectivity. We argue, in consonance with a previous analysiJ^, that 
all phases are characterized by gaps to local excitations. The spin glass states we report arise 
from random long ranged loops which frustrate Neel ordering despite the lack of randomness in the 
coupling strengths. 



I. INTRODUCTION 

The study of quantum antiferromagnets has proven 
among the most enduring themes in modern condensed 
matter physics. The interplay between frustration and 
quantum fluctuations leads such systems to exhibit a va- 
riety of interesting ground states. In this context the 
lattice models constructed by Affleck et. afi^ are par- 
ticularly useful for they build in a great deal of both 
these effects using simple local projectors, which allow 
their ground states to be determined analytically. These 
AKLT models have spins given by S = |M, where M 
is a positive integer, and z the lattice coordination num- 
ber. In principle, they may be defined on any graph, 
but in practice one usually maintains fixed connectivity 
in order to have the same spin on each site. The asso- 
ciated ground states have the added feature that their 
wavefunctions can be written in Jastrow (pair product) 
form, which allows us to view their ground-state proba- 
bility densities as Boltzmann weights corresponding to a 
nearest neighbor Hamiltonian for classical vector spins on 
the same lattice. Using this unusual quantum-classical 
equivalence one can understand many properties of the 
states by studying the associated classical model. 

The initial construction of the AKLT models was mo- 
tivated by the search for quantum disordered states in 
low dimensions. This works only too well: in d = 1 
and d = 2 the mapping to finite temperature classical 
models discussed above ensures, by the Mermin- Wagner 
theorem, that all cases lead to quantum paramagnetic 
or valence bond solid ground states. In d > 2 this is 
no longer true and a computation is needed to decide 
which models order and which do not. In a recent paper, 
Parameswaran, Sondhi and Arovas^ showed via Monte 
Carlo simulations and mean-field arguments that AKLT 
models on the diamond and pyrochlore lattices exhibit 
quantum-disordered ground states for small spin sizes 
while on the cubic lattice all models exhibit Neel order. 

In this paper we take this exploration of higher dimen- 



sional AKLT models in a different direction — we study 
them on locally tree-like lattices of fixed connectivity 
which are known to physicists as Bet he lattices. Here we 
shall consider two physically distinct systems. The first, 
is the Bethe lattice constructed as the limit of a family 
of Cayley trees. This construction yields a system with 
a finite surface to volume ratio and without loops. The 
second is a typical member of the ensemble of random 
graphs of fixed connectivity. These graphs are locally 
tree-like in the thermodynamic limit; however they also 
have long loops of logarithmically divergent size. These 
loops of both even and odd lengths introduce topologi- 
cal frustration into the system. The two constructions of 
locally tree-like lattices yield different physics. 

For the infinite Cayley tree, we exhibit an exact solu- 
tion using the quantum-classical correspondence. Specif- 
ically, we use a generalized transfer matrix technique to 
obtain exact solutions for various statistical quantities in 
the ground state of the tree. We note that the AKLT 
model on the Bethe lattice has been studied before di- 
rectly within the quantum formalisnP; we suspect that 
readers will find our solution simpler. We find one quan- 
tum disordered state (M = 1 on the z = 3 tree) and 
two that are critical (M = 2 on 2; = 3 and M = 1 
on 2: = 4), in that the correlation functions decay ex- 
ponentially at precisely the rate required to balance the 
exponential growth of the graph. All other cases exhibit 
Neel order. We address the question of whether the bulk 
excitations are gapless in cases when the AKLT wave- 
function has critical or Neel correlations. We find, per- 
haps surprisingly, that the system is always gapped to 
local excitations, and that the only gapless excitation is 
a global one connecting the different broken-symmetry 
ground states. We connect this to related work in Ref.[TJ 
in which we conjecture that this is a generic feature of 
symmetry breaking quantum models on the Bethe lat- 
tice, related to the spectrum of the graph Laplacian. 

On random graphs of fixed connectivity, Neel order- 
ing in the companion classical model is frustrated by the 



2 



presence of the long loops. To study this case, we appeal 
to the cavity techniques familiar from the theory of classi- 
cal disordered systems. These have been applied recently 
to a variety of discrete statistical mechanical problems on 
random graphs and there is much evidence that the (ap- 
proximate) techniques are on solid ground. From this 
analysis, we conclude that for 2; < 10 there are disor- 
dered states at small spin and spin glass ground states 
at larger spin as well as a couple of cases where the state 
is critical. For z > 10 all AKLT models have ground 
states with spin glass order. By spin glass order, we 
mean states with fixed but randomly oriented local mag- 
netizations and that the set of such states is larger than 
those connected by global rotations alone. We argue that 
the spectrum of local excitations above the pure states 
in this set is again gapped. 

Of our various results we would especially like to flag 
these last mentioned. The nature of quantum glass 
phases is a subject of much interest - especially as to 
how much of the elaborate framework of the classical sub- 
ject may be lifted into the quantum world. The AKLT 
construction provides a direct line of approach to this 
problem and does so using Hamiltonians without random 
couplings but from graph disorder alone. 

This paper is organized as follows: in Section [llj we in- 
troduce the AKLT model on an arbitrary graph via the 
Schwinger boson formalism. We proceed to construct a 
companion classical model that captures the structure 
of the ground state wavefunction by introducing a ba- 
sis of SU{2) coherent states. In Section III, we develop 



transfer-matrix technology to solve the companion clas- 
sical model on the (bipartite) Bethe lattice exactly, and 
obtain the transition temperature and correlation func- 
tions in the paramagnetic and Neel-ordered phases. In 



Section IV we investigate the energy gap using a varia- 
tional ansatz for the excited states. Finally, in Section [V| 
we consider the extension of this analysis to the spin glass 
transition expected on regular random graphs and con- 
sider some of the quantum consequences of the classical 
glassy phase. 



II. AKLT STATES: A BRIEF REVIEW 

The central idea of the AKLT approach^ is to use quan- 
tum singlets to construct correlated quantum-disordered 
wavefunctions, which are eigenstates of local projection 
operators. One can then produce many-body Hamilto- 
nians using projectors that extinguish the state, thereby 
rendering the parent wavefunction an exact ground state, 
often with a gap to low-lying excitations. A general mem- 
ber of the family of valence bond solid (AKLT) states can 
be written compactly in terms of Schwinger boson^ 



(1) 



(ij) 



per site is given by zM, where z is the lattice coordina- 
tion number, and the resultant spin on each site is given 
hy S = ^zM. Given any regular graph, the above con- 
struction defines a family of AKLT states labeled by the 
size of their spins S = ^zM. For more details regarding 
this construction and the corresponding Hamiltonians, 
see Ref. 4, 

The AKLT states have a convenient representation in 
terms of SU(2) coherent states, as first shown in Ref. [6l 
In terms of the Schwinger bosons, the normalized spin- 
s' coherent state is given by \n) = {p\)~^^'^ {^iJ^\Y I 0)^ 
where p = 25, with z = {u ^ v) CP^ spinor, with u = 
cos(^/2) and v = sin(^/2) e*'^. The unit vector n is given 
by = z^a^z, where a are the Pauli matrices. In the 
coherent state representation, the general AKLT state 
wavefunction is the pair product ^({n^}) = ({n^}|^) = 



Yl{ij){^i^j ~ ^i^j)^' Following Ref. [6] we may write 
|^({ni})p = exp(— /3i^ci) as the Boltzmann weight for a 
classical 0(3) model with Hamiltonian 



(ij) 



1 



(2) 



This assigns M singlet creation operators to each link 
(ij) of an underlying lattice. The total boson occupancy 



at inverse temperature P = M. All equal time quantum 
correlations in the state |^) may then be expressed as 
classical, finite temperature correlations of the Hamilto- 
nian H^^. 

Some immediate consequences of this quantum-to- 
classical equivalence were noted in Ref. |6l On one 
and two-dimensional lattices, the Hohenberg-Mermin- 
Wagner theorem precludes long-ranged order at any fi- 
nite value of the discrete quantum parameter M. In three 
dimensions, there is no a priori reason to rule out long- 
range order. In fact, as shown in Ref. SI the simple cubic 
lattice has no quantum-disordered states at any M, while 
the diamond lattice has a single such state for M = 1; 
on frustrated lattices such as the pyrochlore, such states 
are believed to exist for many values of M. 

As is evident from ([T]) , we may define the AKLT states 
on an arbitrary graph; if the graph has fixed connectivity 
z, then the resulting model has the same spin on each site. 
On graphs with a boundary, this is not automatic, since 
the boundary sites will have fewer neighbors z\ There are 
several ways to deal with this boundary effect. The first 
is to work with a system with a lower spin on the bound- 
ary: in that case = ^z'M. The quantum state of this 
non- homogeneous system is unique. Another option is to 
add {z — z^) additional Schwinger bosons of either flavor 
to the edge sites; there is not a unique way in which to do 
this, leading to a multitude of degenerate ground states 
classified by the state of the boundary spins. When trans- 
lated to the companion classical model, the latter option 
can be viewed as connecting each of the boundary spins 
to {z — z') fixed spins, each with an orientation speci- 
fied by the behavior of an independent spin- 1/2 degree 
of freedom; thus, the different degenerate states of the 
homogeneous AKLT model on a graph with boundary 
can be understood by choosing different fixed boundary 
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conditions for spins in an additional, outer ring of leaf 
spins. Finally, one can opt to get rid of the boundary 
by, for instance, taking periodic boundary conditions on 
a Euclidean lattice. As we will discuss in Sec. |V| when 
generalized to tree-like graphs, this approach leads to a 
spin glass phase in the companion model and thus pro- 
vides a non-trivial new quantum spin glass to the AKLT 
phase diagram. 



III. TRANSFER MATRIX SOLUTION OF THE 
CLASSICAL PROBLEM ON TREES 



Before we proceed, we note that on bipartite graphs 
we can perform a gauge transformation by flipping every 
spin at odd depth to obtain a ferromagnetic model. This 
gives us 



-M ^ log 



1 -\- hi ■ hj 



(4) 



The classical Hamiltonian that describes the properties 
of the ground-state wavefunction of an AKLT model with 
singlet index M, written in the basis of SU{2) coherent 
states is given by 



/3i7ci = -M^log 



1 



(3) 



As usual in the treatment of tree models, we consider 
the statistical state V^^(no) (marginal distribution) of a 
cavity spin no at the root of a branch of the tree. This 
unnormalized distribution can be found in terms of the 
cavity states of its z — 1 neighbors by summation: 



where 



T(no,ni) 



,/31og[i±%:^ 



1 + no • ni 



(6) 



is the transfer matrix of the AKLT model and 

M(no; n[, • • • n^_i) = S{ho - n[) • • • S{ho - n^_i) (7) 



is the merge matrix. 

The merge matrix M defines a multilinear map from 
the z — 1 state spaces of the neighbor spins to the state 
space of the root. This lifts naturally to the appropriate 
complexified tensor product spaces^ and thus we will 
find it natural to write the merge and transfer operations 
abstractly using Dirac notation: 



M 



I 



DhoDh[ ' ••Dn'^_^ S{ho - h[) ■ - ■ S{ho - K-i)\^o){n[ \ ■ ■ 
Dh |n)(n| • • • (n| 



(8) 



and, 

T = j DnDn'T(h,n')\fi){fi'\. (9) 

Thus, equation ([5| becomes 

IV^O) = M (TlV^i) . • • T\r-^)) . (10) 

We now focus on the stability of the paramagnetic state 
against Neel ordering. Hence, we have to use boundary 



conditions that are consistent with this kind of ordering. 
As discussed at the end of Section [Hj one can either use 
free boundary spins with lower or connect the bound- 
ary spins to some fixed additional spins: in the latter 
case, the additional spins must all have the same orien- 
tation to allow the Neel ordering. This remark is particu- 
larly important because, as we will see later in Section [V} 
on a random regular graph the boundary conditions on 
any given tree-like subregion are fixed self-consistently. 
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and in general are not consistent with Neel ordering, lead- 
ing to a disordered spin glass state. Assuming uniform 
boundary conditions, we obtain the unique state at depth 
d — 1 hy merging the z — 1 states at level d using the T 
and M operators: 



\d-l) =M{T\d))®^'-^\ 



(11) 



The natural basis to work in is that of states with 
definite angular momentum, i.e. states |/m), which are 

I 



eigenstates of the angular momentum operators 1/^,1/^. 
In the coordinate basis, these are simply the spherical 
harmonics, and as shown in Appendix |Aj they are eigen- 
states of the transfer matrix with eigenvalue A^. 



It remains for us to understand exactly how the merge 
operation acts in the angular momentum basis. If we in- 
sert resolutions of the identity in the angular momentum 
basis into (|8|, we obtain 



M = 



J Dh\h){h\---{h\ 

lo,mQli,mi Iz — i^fTLz — i 



,mo){lo,mo\n) {{n\h,mi){h,mi\ (g) • • • (g) {h\lz-i,mz-i){lz-i,mz-i\) 



E E--- E l«o,mo)((;i,mi|®...®(?,_i,m,_i|)|Dnl^7*(n)y™^(n)...l^^:^^ (12) 

lo,moli,mi lz-i,mz-i 



r 



For the case z = 3, the integral in (12) is simply the 



Clebsch-Gordan coefficient that characterizes the fusion 
of two SU(2) spins. For higher values of this is the 
appropriate generalization of the Clebsch-Gordan coeffi- 
cient describing the fusion of {z — 1) SU{2) spins. Thus 
we see that the merge operation, when written in the an- 
gular momentum basis, has a natural interpretation as 
the fusion rules for the 0(3) symmetry group. 

The paramagnetic state - here represented in Fourier 
space by the |00) state - is always a fixed point: it is 
an eigenstate of the T-matrix, and the fusion of any 
number of |00) states is again a |00) state. We pro- 
ceed via linear stability analysis: we introduce a per- 
turbation into a state that is not uniformly weighted on 
the sphere, and see if this grows or shrinks under the 
iteration procedure. We note that we can decompose 
any such state into spherical harmonics, and so we write 
\d) = 1 00) ^ e^'^Q^^cim\lm) into ([iT]) to obtain 



1^-1) = 



M I Ao|00)+e Xicim\lm) 



where we have used the fact that fusing any number of 
1 00) states with an |/ m) state results in an |/ m) state. We 
renormalize to leading order and find that the iterated 
state is, to linear order 



\d-l) = |00) + e(z-l) 



E 



ci^\lm) + Oie^) (14) 



The perturbation is irrelevant (shrinks under iteration) 
if the coefficient of the linear term is less than 1, and rel- 
evant if it is greater than 1. The critical point is reached 
when 



hi 

Ao 



1 



(15) 



for any Using the temperature dependent expression 
(A4) for the A/, one can show that the dipole instability 
(r= 1) is the first one encountered as the temperature is 
lowered, and therefore sets the transition temperature. 

Using the results of Appendix |A] (replacing (3 by the 
singlet parameter M), we obtain 



Mr = 



(16) 



We see that for z = 2,3,4, Mc = oo,2,l, while for all 
other values, Mc < 1. Since M must be a positive in- 
teger, we see that for the chain {z = 2) all values of 
M correspond to quantum-disordered states (which fol- 
lows from the Mermin- Wagner theorem and the original 
AKLT result^) whereas for z = 3, the M = 1 state is 
disordered while the M = 2 state is critical, and finally 
for z = 4, the M = 1 state is critical. Bethe lattices 
of higher connectivity will always have ordered AKLT 
ground states for any value of M. See Fig. [l] 

Finally, we consider the correlation function (no • ^d) 
within the paramagnetic phase. This is given by consid- 
ering the response of (no) to a field on fid - which is the 
same as asking how the dipole / = 1 perturbation propa- 
gates along a chain of length d in a background of trivial 
/ = cavity states. This immediately implies 



(no • fid) (X 



M 



M + 2 



(17) 
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Notice that this imphes that the naive correlation length 
never diverges - as usual with tree models, phase transi- 
tions occur when the correlation decays at the same rate 
as the growth of the lattice. For a slightly more detailed 
calculation, see Appendix [B) 

As an aside, we note that we can use the same gener- 
alized transfer matrix technique to obtain the transition 
temperature for the Heisenberg model, a result first ob- 
tained by Fishei!^ using a different method. This serves 
as a test of the technique proposed here. 



IV. VARIATIONAL BOUNDS ON THE GAP 

We now perform a variational computation of the gap 
to excitations in the critical model, similar to the Single- 
Mode Approximation (SMA) discussed in Ref. 6, The 
central idea of the SMA is to construct an excitation 
orthogonal to the ground state by acting on it with a local 
operator, and then to reduce the energy of this excitation 
by delocalizing it, thereby decreasing its kinetic energy. 
A variational bound on the energy gap is given by 



< A < Asm A 



{^sma\H — Eq\^sma) 



smaV^sma 



(18) 



This approach is designed to optimize the energy due 
to the off-diagonal matrix elements in the excited sec- 
tor, which will be proportional to the usual graph Lapla- 
cian for a nearest neighbor model. On the Bethe lattice, 
the spectrum of the Laplacian is unusual. As argued in 
Ref. 1 , there is necessarily a gap to hopping excitations on 
tree-like graphs despite the existence of symmetry related 
ground states. Thus, in some sense the SMA calculation 
is doomed to failure as it will never be able to close this 
gap. Nonetheless, it is interesting to see how this plays 
out in an exact model. 

We consider a rooted Cayley tree, and in order to re- 
strict ourselves to studying excitations confined to the 
bulk as the size of the tree grows we suppress excitations 
far from the center using an infrared regulator A. We 
therefore study variational wavefunctions of the form 



N 



(19) 



where ux{i) is a function only of the depth pi of site i 
referenced to the root of the Cayley tree. The SMA gap 



^SMA 



lim 



lim 

N^oo 



{X\H-Eo\X) 



(A|A) 



(20) 



After some algebra, we may write (with the understand- 
ing that we always take N to infinity before taking A to 
zero): 



^SMA = lim 

N^oo 



(21) 



Asserting homogeneity of the graph, and noting that 
the Hamiltonian is a function only of Si • Sj where i and 
j are nearest neighbors, we may re-express this as 



^SMA = lim 

N^oo 
A^O 



fx 



(22) 



where / = ^([S'f - [H,Sf - and we denote 

by A the adjacency matrix of the Cayley tree; the corre- 
lations in the ground state are always exponential, and 
go as 7'*"-^', where 7 -^-^ from below as we approach 
criticality. We have therefore reduced the problem of 
the SMA on the Bethe lattice to understanding (i) the 
spectrum of the graph Laplacian (the adjacency matrix 
up to a sign) and (ii) the behavior of the ground-state 
correlations. 

Our choice for a variational ansatz is to take u\{i) = 
e~^Pi . This is motivated by the fact that the number 
of sites at a given distance from the center grows expo- 
nentially, and therefore in order to remain in the bulk of 
the tree, we need to cut off the wavefunction faster than 
exponentially.^^ We perform the summations by convert- 
ing the sum over sites into a sum over depths, approxi- 
mating the sums by integrals and using steepest-descent. 
We find that, at criticality, the gap is nonvanishing: 



\G f ^ 

^SMA ^ 



2^fz 



1 



1 



5(^-1) 



(23) 



Excitations constrained to live in the bulk are there- 
fore always gapped, even at criticality. The factor of 
z — 2^/ z — 1 is precisely the spectral gap for bulk excita- 
tions on a Cayley tree^. A state where ux{i) is indepen- 
dent of position must be gapless in the broken-symmetry 
phase of the model, since it connects the different broken- 
symmetry ground states. We cannot recover this state 
through a correctly regulated calculation in the chosen 
order of limits, however. 

There is a straightforward physical argument for this 
gap. By our choice of variational ansatz, we cut off the 
excitation at some depth D ~ 1/a/A. On a Euclidean 
lattice of dimension (i, this costs a surface energy ex D^~^ 
which is normalized by the weight of the wavefunction 
(x D"^. As i:) ^ 00 (A ^ 0), the SMA gap therefore 
vanishes as . On tree-like lattices, both the surface 
area and the bulk normalization scale as {z — 1)^; the 
boundary is always a finite fraction of the bulk. Thus, 
the ratio remains finite as I) ^ 00 and the gap cannot 
clos^. 



V. AKLT MODEL ON REGULAR RANDOM 
GRAPHS: FRUSTRATION AND THE SPIN 
GLASS STATE 

We now consider the same model on a regular random 
grapJ^ of connectivity z. The ensemble of these graphs is 
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constructed by assigning uniform probability to all pos- 
sible graphs of vertices, such that each vertex is con- 
nected to exactly z links. There are several reasons why 
statistical models defined on this ensemble of graphs are 
interesting: 

1. A central property of this ensemble^ is that typ- 
ical lattices are locally tree-like; their loops have 
a length diverging logarithmically with the size N 
of the system; this implies that one can develop a 
method to solve statistical models on these graphs 
based on the same recurrence equations that are ex- 
act on trees. This is known as the cavity method^^. 

2. Despite being locally tree-like, they do not have any 
boundary, all sites playing statistically the same 
role (in the same way as periodic boundary condi- 
tions impose translation invariance on a finite cubic 
lattice). Moreover, the free-energy of regular ran- 
dom graph models is self- averaging with respect 
to their random character in the thermodynamic 
limit. In other words for large enough N a sin- 
gle sample is a good representative of the ensemble 
average. 

3. Typical graphs are characterized by many large 
loops of even and odd length; this strongly frus- 
trates the antiferromagnetic ordering, which gives 
way to a spin glass phase instead^^. 

In the following we will be particularly interested in the 
implications of the last point for the quantum problem. 

The reasoning outlined in section [TT| clearly applies to 
the random graph model, whose AKLT ground state is 
therefore described by a classical Hamiltonian of the form 
(|3|, where the pairs (z, j) are connected by a link of the 
random graph. The main difference between the tree 
model and the random graph model is that the recur- 
rence equation ([5| now does not hold for the full graph; 
it only holds for a tree-like subregion of the graph, and 
has to be initialized using the boundary values of the 
ip'^{hi) that are determined by the summation over the 
rest of the graph. In other words, the recurrence on 
the subregion is initialized from random self-consistent 
boundary conditions, determined by the rest of the sys- 
tem: these boundary conditions are not consistent with 
Neel ordering, which is therefore frustrated, as discussed 
in Section [TTT| above. However, since the tree-like subre- 
gions grow in size when TV ^ oo, equation ([5| must be 
iterated a very large number of times. One can classify 
the different phases of the system by studying its fixed 
point^i^ 

To calculate the stability of the paramagnetic solution 
against spin glass ordering we observe that the spin glass 
transition is signaled, as usual, by the divergence of the 
classical spin glass susceptibilit}'^^: 



XSG 
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FIG. 1: Phase diagrams for AKLT model with singlet pa- 
rameter M on tree- like lattices with coordination z. On the 
Cayley tree the transition is from paramagnetic (PM) to Neel- 
ordered (AF) phase at the solid blue line with no spin glass. 
On regular random graphs the transition is from the param- 
agnetic to spin glass ordered (SG) phase at the dashed green 
line — there is no antiferromagnet. The models with Bethe 
lattice critical correlations are labeled with large dots. 



The details are discussed in Appendix |B] the result is 
that xsG is finite if for all / 



Az 1 



Ao — 1 



(25) 



Once again the instability originates in the / = 1 sector 
and occurs at 



Msg 



(26) 



(24) 



We see that Msg = oo for z = 2, Msg = 4.828 for z = 3 
(hence the system is a spin glass for M > 5), Msg = 2 
for z = 5 and Msg = 1 for z = 10 (in these cases the 
system is critical), and it is smaller than 1 for any z > 10. 
See Fig. (T) 

With the Neel-ordered phase suppressed, the quantum 
paramagnet extends further in the z — M plane than 
on the Cayley tree models. Nonetheless, it is clear that 
the paramagnetic solution develops an instability to spin 
glass ordering at large M. What are the properties of 
this low temperature phase? Most of the detailed work 
on classical spin glasses has focussed on discrete models. 
The AKLT mapping provides a classical vector model 
with weakly divergent interactions whose glass phase has 
not yet been studied. Nonetheless, most of the qual- 
itative features of the classical multiple-valley picture 
should still hold and these provide an intriguing scenario 
for the quantum system. The classical Gibbs measure 
decomposes into a collection of clustering pure states 
a = 1...JV with essentially disjoint support. In each 
of these states, the i^l,{ni) - and thus the local magne- 
tizations - are macroscopically different. This strongly 
suggests that the quantum AKLT ground state |^) it- 
self is a superposition over a collection of macroscopically 
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distinct degenerate ground states I^q,) each of which cor- 
responds to one of the classical clustering states. 

While we believe that the above picture holds in gen- 
eral, a rigorous derivation is problematic. In the following 
we will attempt to justify it in more detail and point out 
some of the subtleties that must be dealt with. First, 
there does not yet exist a detailed study of the classi- 
cal vector model^^ - for the purposes of this paper, we 
shall assume that, modulo the 0(3) global symmetry, the 
qualitative behavior is that of the better studied Ising an- 
tiferromagnet on a regular random graph^^. 

By analogy to this model, two different phases exist: 
at high temperatures, the stable phase is a paramagnet 
where '^*(ni) = tp{h) = 1 is the same for all sites, and 
it is the unique fixed point of Eq. ([5|. At low tempera- 
tures, the stable phase is a spin glass, characterized by 
the existence of many non- symmetry related pure states 
labeled by a = 1,...,A/' within which connected spa- 
tial correlators vanish. This is related to the existence 
of many different fixed points of Eq. ([5|, and reflects 
the decomposition of the thermodynamic Gibbs measure 
P = ex.-p{—PHci)/Z as follows: 



P({n}) 



a=l 



dgPa{{g ■ n}) 



(27) 



where the Pa are representative classical pure state mea- 
sures. By integrating over the 0(3) of global rotations 
of spin space, we account for the continuous family of 
symmetry related pure states associated to each repre- 
sentative state a. One can access these representative 
pure states by adding a uniform infinitesimal field to the 
classical model, but as the quantum AKLT state is a sin- 
glet, we prefer to work without explicitly breaking this 
symmetr}ff^. 

The Ising spin glass models with two-body interac- 
tions which have been studied, such as the Sherrington- 
Kirkpatrick modef^ and the random graph antifer- 
romagnetP^, are characterized by a continuous spin 
glass transitiorPI' with a finite collection of pure states 
throughout the spin glass phase^. We shall therefore as- 
sume that this is true of our collection of representative 
pure states. Indeed, all we will need is that A/" grows at 
most polynomially in TV as the thermodynamic limit is 
taken. 



The pure state decomposition (27) has striking conse- 



quences for the structure of the low-energy states of the 
quantum AKLT Hamiltonian. To wit, we argue that 



i*) = E 



a=l 



dgg] 



(28) 



up to exponentially small corrections in the thermo- 
dynamic limit, where the |^q,) can be interpreted as 
a collection of symmetry breaking degenerate quantum 
ground states whose correlations correspond to the clas- 
sical pure states P^. 



We argue this in three parts. First, we assume the 
existence of a collection of quantum states |^q,) such that 



|({nj|vl/„)|2=p„({nj) 



(29) 



and show that the quantum state Eq. ( |28| ) reproduces the 
observables of the classical decomposition of Eq. (27). 
Second, we will show that up to exponentially small 



corrections each of the |^q,) are themselves orthogonal 
ground states. Finally, we address the issue of the exis- 
tence of such states given the classical distributions P^- 
The first part follows the argument of Ref. 15 but we 
rephrase it in terms of density matrices. Consider the 
density matrix of the proposed state Eq. (28): 



p=\^){<i/\ = Y^,/^U^ J dg' j dgg'\^a){^0\g^ 

a, (3 

(30) 

Given any local operator O depending only on spinJ^, 
its expectation value in state l^f) is given by 

(O) = TrpO 

= J2V^W J dg' f dgTv [5'|*«)(*^|5ta 

a, (3 

= EV^W / dg' J dg{^^\g^dg'\^^) (31) 

a, (3 

We now argue that the interference term is negligible. 
That is, 

(^^|^ta^/|^^) = Sa^5gg.{^Ila\g^dg\^a) (32) 

up to exponentially small corrections in the thermody- 
namic limit. This follows from the observation that l^^a) 
and 1^/3) have macroscopically distinct magnetization 
patterns that completely break the 0(3) symmetry. In 
particular, the classical configurations {n^} on which the 
wavefunction |^q,) is concentrated have extremely small 
weight in any other wavefunction P ^ a - this remains 
true even with arbitrary global rotations allowed between 
them. lia = 13 but g and g' differ, then the configurations 
with weight are again macroscopically distinct by virtue 
of the net global rotation g~^g' between them. The fact 
that the observables are local and have bounded matrix 
elements does not modify these assertions. Since M is fi- 
nite for our antiferromagnetic model, the finite sum over 
the exponentially small corrections does not modify the 
result: 

{6) = Y.Wa j dg{^aWdg\^a) (33) 

Inserting a complete set of coherent states and using 
Eq. ([29| reproduces the classical distribution Eq. (27). 
The probability of finding the quantum system in a state 
a (with respect to the state |^)) is the same as that of 
the classical problem: both are given by Wa- 
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Furthermore, the l^^a) must have exponentiahy smah 
energy with respect to the quantum Hamiltonian. Since 
(H) = in the AKLT state, using equation (33) and the 
rotational invariance of H, we have that 



0={H) 



(34) 



up to exponentiahy smah corrections. Since each of the 
terms in the sum is nonnegative, it fohows that 



(vl>„|i/|vl>„) < 0(e-^) 



(35) 



Thus, the |^q,) are a cohection of nearly orthogonal, 
nearly zero energy states, each of which generates a fur- 
ther continuous collection of such degenerate states under 
the action of 0(3). 

Finaly, we turn to the slightly thorny question of 
whether states satisfying Eq. (29) exist. The problem 



is that the coherent state basis is overcomplete for any 
fixed spin size S = zM/2 and we cannot necessarily find 
quantum states which have given expansions in this ba- 
sis. That is, a priori we cannot simply set {{hi}\'^a) = 

Pa{{ni}) and know we have a well-defined quantum 
state for spins of size S. In the large spin limit, there 
is no problem as the coherent states become a complete, 
rather than overcomplete, basis. This coincides with the 
zero temperature limit of the classical companion model 
and thus the pure states I^q,) may be identified with the 
(many degenerate) minima of the energy function (|3|. At 
finite M, we cannot find such finely localized states in the 
coherent state representation - the most localized state 
has solid angular scale ~ 1/M - but the finite temper- 
ature fiuctuations around the classical minima will also 
smear the Pa at a similar scale. 

On the other hand, we already used above that the 
states have disjoint support, up to exponentially small 
corrections in N. For a given classical configuration 
{n^}, only one state contributes to ({ni}|^) significantly. 
Thus, in a given region of classical configuration space, 
({ni}|^) coincides with one of the {{hi}\'^a)^ and con- 
versely, each of the {{ni}\'^a) can be seen as the restric- 
tion of the full ({n^}|^) to that region. Since the 1/M 
smoothing is local in configuration space, it is safe to 
assume that the states ({n^H^Q,) are as smooth as the 
original ({ni}|^), therefore |^q,) can be defined without 
ambiguity. Based on the above arguments, we think it 
likely that at least an approximate pure state decompo- 
sition of the form proposed above can be found even at 
finite M. 

In summary, we obtain the following picture for the 
low energy spectrum in the spin glass phase: the non- 
clustering paramagnetic AKLT ground state |^) can 
be decomposed in a superposition of several almost de- 
generate states, whose energies are of order exp(— A/"). 
These states enjoy the clustering property (vanishing of 
connected correlations) and are characterized by amor- 
phous order (the local magnetizations are different in 
each state). It would be nice to check this scenario ex- 
plicitly by means of exact diagonalization. 



It would be very interesting to obtain more detailed 
information on the spectrum of such spin glass Hamil- 
tonians. For instance, a natural question is whether 
there is an energy gap between the degenerate low-lying 
spin glass states and the excited states. Indeed, we 
expect Goldstone (or Halperin-Saslow) modes^^ associ- 
ated with twisting of the amorphously magnetized states 
l^o;). While it is difficult to explicitly construct a coarse- 
graining procedure to produce an effective theory of such 
modes on a tree-like graph, one usually expects that such 
a theory applies in sufficiently high dimensions'^. Insofar 
as the effective theory is an elastic hydrodynamics living 
on a tree-likegraph, the corresponding modes should re- 
main gapped^. This suggests that the low energy spec- 
trum is indeed gapped in any given pure state sector. 



VI. CONCLUDING REMARKS 

In this paper, we extend the study of AKLT models to 
locally tree-like graphs of fixed connectivity by exploiting 
the quantum-classical mapping of the associated wave- 
functions. On the infinite Cay ley tree, we recover the 
results obtained in Ref. 5 . We find that the Bethe lat- 
tice possesses the peculiar property that it is possible to 
choose parameters (for 2; = 3, 4) so that the correspond- 
ing AKLT state is critical. A variational calculation of 
the gap is unable to produce gaplessness, which is con- 
sistent with the arguments of Ref. [T] that this is a general 
feature of locally tree-like graphs: essentially, one cannot 
deform a uniform excitation into long-wavelength rota- 
tions of the order parameter, without jumping a gap in 
the Laplacian spectrum. 

Turning to regular random graphs, we find that the 
companion classical model is unstable to spin glass or- 
dering within a cavity analysis. This is a general feature 
of classical antiferromagnetic models on such graphs, but 
has striking consequences given that the pecularities of 
such mean-field-like glasses should directly transfer to 
the quantum ground state of the AKLT model. This 
provides an alternative route to the study of quantum 
glassy order in tree-like models (see for example Ref. 18]- 
23). We argue that there are now many (nearly) degener- 
ate quantum ground states with macroscopically distinct 
magnetization patterns, but that there remains a gap to 
Halperin-Saslow waves for geometric reasons analogous 
to the simpler case of the antiferromagnet. 

There are several avenues for future research. One ob- 
vious direction is to study the classical vector spin glass 
and the corresponding classical measure. In a different 
vein, we observe that the AKLT construction applies at 
a very special point in the space of quantum Hamilto- 
nians. To what extent do the features of the quantum 
AKLT glass extend to regions proximate to this exactly 
solvable point? Ideally, the AKLT glass would capture 
the essential features of a broader range of quantum spin 
glasses, playing a role reminiscent of that played by the 
S = 1 AKLT chain in relation to the Haldane phase. 
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(Bl) 



can be rewritten, by taking the average over the random 
graphs and using translational invariance, as 



XSG = ^J^d[{no • nd)f 



(B2) 



d=0 



where JVd is the number of sites at distance d from a 
reference site. The sum is convergent as long as 



Appendix A: Transfer Matrix for the AKLT Model 



lim (A/-d)i/'^[(no • fid)]'/^ < 1 (B3) 



The transfer matrix may be thought of as a map be- 
tween functions defined on the sphere: 



(r/)(n) 



dn' 
47r 



T{h,h')f{n') 



(Al) 



where n, G S^-\ In our case, as in most cases of 
interest, the transfer matrix is a rotational scalar and we 
can work in the angular momentum basis. The eigenvalue 
must depend only on the eigenvalue / and not on the 
Lz eigenvalue m. It therefore suffices to solve the problem 
in the case m = 0. 

We wish to solve the eigenvalue equation: 



... f Dh' 



T{n,n')fi{n') 



(A2) 



Since the kernel T depends only on n • = cos^, we 
work in polar coordinates with the z-axis along ft and 
substitute x = cos 9 to obtain 



) = \j\xT{x)fi{x) (A3) 



A natural guess for the eigenfunctions is that they are 
Legendre polynomials. The transfer matrix in our case 
is T{x) = (^^)^. Using standard identities. 



A, 



^,AKLT 



[r(/3 + i) 



■l + l)V{f} + l + 2) 



(A4) 



A similar discussion for the Heisenberg model for arbi- 
trary N and the case of SU {N) and Sp{N) groups, may 
be found in Ref. EH 



Appendix B: Stability against spin glass ordering on 
a regular random graph 

We follow closely the derivation of Ref. ^25^, Appendix 
A. In the thermodynamic limit, the spin glass suscepti- 



Note that (A^^) 



l/d 



2; — 1 for large d. In the param- 



agnetic phase, (no • ^d) is given by the response of (no) 
(the root) to a small magnetic field coupled to n^^, a leaf 
at distance of a tree whose other nodes are in the 
paramagnetic state |00). Hence we get (repeated indices 
summed): 



(no • fid) cx 



dhl, 



I 



dhi, 



(B4) 



Clearly the term that gives the exponential dependence 
on d is the variation of '^^(no) with respect to hd- Using 
the recursion relation ^ we can rewrite it as 



#'(no) 
dhd 



dhi ■ • • dhd 



#0(no) #^(ni) di/j'^ihd) 
dil)^{hi) dil)'^{h2) 



dhd 



and the exponential dependence is related to the eigen- 



value of the transfer matrix 



= T{hd,hd^i). 



di;d+^hd+i) 

These can be obtained by repeating the analysis of sec- 
tion |TTI1 Indeed if we use the same ket notation and write 
the variation S\d -\- 1) = e X][^o m ^im\^ ^) 
results of section [Till to obtain 



A. 



(B5) 



The absence of the factor {z — 1) with respect to Eq. (14) 



is due to the fact that here we are only varying one of 
the neighbors of a given spin, the neighbor on the path 
linking the root to the given leaf at distance d. The 
relevant eigenvalues are therefore Xi/Xo and we obtain 
the condition 



(z-l)max{(A^/Ao)'}<l 
for the convergence of XSG- 



(B6) 
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A few comments on the nature of the classical statistical 
state space of a vector spin are in order, as we have so 
cavalierly complexified and tensored it into a much more 
quantum mechanical looking system. Physical cavity dis- 
tributions '0(n) must be real, normalizable, nonnegative 
functions on the sphere. By its construction as a marginal- 
ization (summing out) procedure, equation |To| must pro- 
duce such a physical output given physical inputs, even 
though we have extended it over C. A more important sub- 
tlety arises in the normalization of states - the standard LF' 
norm associated with the Dirac inner product is not neces- 
sarily 1 for a properly normalized probability distribution. 
Since the probabilistic norm is incompatible with the 
Hilbert space structure, it is much simpler to work with un- 
normalized vectors and keep in mind that a probabilistic 
interpretation only applies in the standard basis. 
In principle, an exponential regulator with a sufficiently 
fast decay constant also works. The calculation proceeds in 
a similar fashion but the interpretation is more complicated 
and no more enlightening. 

See Ref. 10 for a general discussion and Ref. [13| for the 
explicit computation of the phase diagram of a classical 
Ising antiferromagnet. 

Although there is a replica symmetric treatment of a re- 
lated model with the additional complication of random 
interactions. See Ref. [26l 

Other mean field spin glass models, such as the p-spin 
model, can have a number of states scaling exponentially in 
system size. These models show a discontinuous spin glass 
transition. However, antiferromagnetic models with two- 
body interactions usually do not show this phenomenol- 
ogy, therefore we will not investigate this transition in this 
paper. 

See Ref. 14, and in particular the reprint on page 226, for 
a more detailed discussion of this delicate statement. 
Such observables are diagonal in the coherent state basis 
and therefore their correlations follow from the classical 



